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Phylogenetic Relationships of the Genus Meretrix (Mollusca: 
Veneridae) Based on Mitochondrial COI Gene Sequences 
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Abstract: Fifteen sequences from the mitochondrial cytochrome c oxidase subunit I gene (COI) were determined for 
4 species of the genus Meretrix, with the homologous sequences of M. petechialis obtained from the GenBank data library. 
The alignment length of the sequences was 574bp after excluding ambiguous sites, including 93 parsimony informative 
sites. In the fragments, the percentages of A, T, C and G were 21.15%, 44.71%, 14.05% and 20.09% respectively. There 
were 12 haplotypes identified: 4 M. meretrix, 2 M. lamarckii, 3 M. lusoria, 1 M. lyrata and 2 M. petechialis. Furthermore, 
it was revealed that M. meretrix, M. petechialis and M. lusoria shared some haplotypes. Phylogeny trees were 
reconstructed by Maximum-parsimony (MP) and Bayesian method using Cylina sinensis as the outgroup. Our results 
indicated that M. lusoria, M. petechialis and M. meretrix are closely related species. This is in accordance with the 
viewpoint that M. petechialis and M. lusoria should be treated as a junior synonym of M. meretrix. 
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Members of the genus Meretrix are highly among systematists in the world. 
appreciated table food for internal and export markets, Traditionally hard clams are mainly identified based 
fetching a high price. Some species have already become on visible shell characters. Commonly, taxonomic 
regional economic mainstay in coastal areas of China. ambiguities exist due to morphological variability. 
However, previous research on species of the Meretrix Modern taxonomic work includes analysis of a host of 
mainly focused on aquaculture and fishery and paid little other traits, including anatomy, physiology, behavior, and 
attention on its classification. Thus, the taxon of the genetics. Citations of clam studies include many 
Meretrix is still in argument over several discrepancies synonyms, indicating ambiguities in Meretrix species 
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identification. Accurate species identification is very 
important in the case of multiple morphology-based 
classification and determining similar species for 
aquaculture management, biodiversity studies, and 
population dynamics. 

The mitochondrial cytochrome c oxidase subunit I 
gene (COI) shows distinct divergence and provides 
valuable information in species identification to 
complete taxonomic data and validation of systemic 
position, phylogeny (Machordom et al, 2003; Smith et al, 


2004; Donald et al, 2005). In the present study, the COI 
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gene is selected to elucidate taxonomy and phylogenetic 
relationships among 5 Meretrix species from the Chinese 
coast. 


1 Materials and Methods 


1.1 Samples and DNA extraction 

Geographical origins of the Meretrix used here are 
shown in Tab. 1. All the clams were identified 
morphometrically, and then fifteen of each species’ 
adductor muscle tissue samples were collected and 


preserved in 100% alcohol. 


Tab. 1 Collection locality and GenBank accession numbers of the specimens studied 





Genus Species and haplotype number identifying the sequence Locality GenBank accession numbers 

Meretrix M. meretrix 1 Jiangsu, China FJ434675 
M. meretrix 2 Jiangsu, China FJ434676 
M. meretrix 3 Jiangsu, China FJ434677 
M. meretrix 4 Jiangsu, China FJ434678 
M. lusoria 1 Jiangsu, China FJ434679 
M. lusoria 2 Jiangsu, China FJ434680 
M. lusoria 3 Jiangsu, China FJ434681 
M. lyrata Hainan, China FJ434682 
M. lamarckii 1 Hainan, China FJ434683 
M. lamarckii 2 Hainan, China FJ434684 
M. petechialis 1 AB280785* 
M. petechialis 2 AY874530* 

Cylina C. sinensis 1 Jiangsu, China FJ434685 
C. sinensis 2 Jiangsu, China FJ434686 


* Sequences downloaded from GenBank. 


Total genomic DNA was extracted from the stored 
tissue samples by the standard Proteinase-K/Phenol- 
Chloroform-ethanol method. Before incubation the 
samples were soaked in ultra pure water for 2 days. 
Scissored tissue samples were re-suspended in 400 uL 
digestive system, which contained Tris (pH 8.0) 0.01 
mol/L, EDTA (pH 8.0) 0.1 mol/L, NaCl 0.05 mol/L, SDS 
1% and 10pL Proteinase K, and was incubated at 52°C 
for 12h. The digested samples were phenol-extracted, 
ethanol-precipitated once more, and redissolved in 10 
mmol/L Tris-HCl, pH8.0. The concentration of DNA was 
estimated using a UV spectrophotometer. The DNA was 
diluted to a final concentration of 50ng/uL. To avoid 
DNA damage caused by multigelation, DNA samples 
were stored at 4°C. 

1.2 PCR amplification 

The partial region of the mitochondrial COI gene 
was amplified from 15 individuals of each species using 
the universal invertebrate primers COI-F 5'-GGTCAA- 
CAAATCATAAAGATATTGG-3' and COI-R 5'-TAA- 
ACTTCAGGGTGACCAAAAAATCA-3' (Folmer et al, 
1994) synthesized by Shanghai Sangon Company, China. 


The polymerase chain reaction (PCR) is employed 
in this study to amplify the COI region. A total 30uL of 
each PCR reaction comprised the following: PCR buffer 
(10X) 3.0pnL, MgCl, (25 mmol/L) 2.0 uL, dNTPs (2 
mmol/L) 2.0 uL, forward primer (10 umol/L) 1.0 uL, 
reverse primer (10umol/L) 1.0uL, Taq polymerase 0.2 
uL, DNA (50ng/uL) 1.0uL, ultra pure water 19.8 uL. 
The Hot-Start PCR was 
denaturation of 5min at 94°C followed by 30 cycles of 
denaturation for 30s at 94°C, annealing at 52°C for 30s 
and an extension of 72°C for 45s. After the completion 


employed with initial 


of 30 cycles, a final extension step of 7min at 72°C was 
performed. The PCR product was then kept at 15°C until 
removed from the machine. The amplified product was 
tested in 2% agarose gel and visualized using the Gel 
Doc system. Those products with distinct and intense 
bands were selected after a series of sequencing 
protocols. The sequencing was done both in the forward 
and reverse directions by Shanghai Sangon Company. 
1.3 Data analysis 

Fifteen COI gene sequences of 4 species in the 
Meretrix (except sequences of M. petechialis) were 
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sequenced, and COI sequences of M. petechialis were 
downloaded from GenBank. The raw nucleotide 
sequences obtained were assembled using SeqMan 
(DNAstar package) and the sequence information 
checked entirely for consistency from both directions 
manually. We examined 574 base pairs (bp) and 
identified 12 haplotypes: 4 of M. meretrix, 2 of M. 
petechialis, 3 of M. lusoria, 1 of M. lyrata and 2 of M. 
lamarckii (Tab. 1). In which, the sequence of M. lusoria 
3 was identical to the sequence of M. petechialis 2, and 
the sequence of M. petechialis 1 was identical to the 
sequence of M. meretrix 2. The collated sequences of 12 
haplotypes were aligned using Clustal X1.8 (Thompson 
on default. No 
insertions/deletions and stop codons were examined. 


et al, 1997) with parameters 
These sequences were translated according to the 
invertebrate mitochondrial genetic code to the expected 
191 amino acids. Nucleotide variation and substitution 
patterns were examined using the software MEGA3 
(Kumar et al, 2004) on the basis of uncorrected pairwise 
genetic distance (p-distance). 

Maximum-parsimony (MP) analyses were 
performed on PAUP 4.0b10 (Swofford, 2002), using 
Cylina sinensis as an outgroup species. Heuristic MP 
Search was performed using the simple addition 
sequence and the tree bisection-reconnection (TBR) 
branch swapping algorithm. Both weighted and 
unweighted analyses were performed. Weighting was 
carried out to equalize contributions from each codon 


position: positions weighted 3: 11: 1 for first, second, 
and third positions, respectively. Levels of support for 
individual relationships were estimated through 1,000 
bootstrap replicates. 

Bayesian inference was done using MrBayes 3.0b4 
(Huelsenbeck & Ronquist, 2001) with the best-fitting 
model TrN+G estimated by Modeltest 3.7 (Posada, 1998). 
Trees saved below the burn-in generations were 
discarded, and a majority-rule consensus of the remains 
were calculated in MrBayes 3.0b4, providing posterior 
probabilities for clades. The MrBayes 3.0b4 was run with 
the following specifications: TrN model with a gamma 
distribution (invgamma), Markov’s chains started from a 
random tree for 400 000 generations, and sampling of the 
Markov chains at intervals of 100 generations. Four 
chains were run simultaneously with the initial 200 
cycles discarded as burn-in. 


2 Results 


2.1 Description of data 

After aligning, using Clustal X1.8, the length of 
COI was 574bp. The obtained COI gene sequences 
contained 149 variable sites and 93 parsimony- 
informative sites. The average base composition were 
21.15% A, 44.71% T, 14.05% C, and 20.09% G. The A 
+ T contents were higher than those of G+C, which is a 
pattern that has been repeatedly seen in the mtDNA of 
Molluscs. 


Tab.2 Nucleotide substitutions and amino acid variability in COI fragments for 12 haplotypes 





All sites Ist codon position 2nd codon position 3rd codon position Amino acds 
Total 574 192 191 191 191 
Invariant 425 158 184 83 150 
Variable 149 34 7 108 41 
Parsimony 93 20 6 67 29 


Nucleotide variation and substitution patterns were 
examined using the software package MEGA 3.0. The 
average values of intraspecific pair-wise sequences 
4.5%. The highest 
divergences were among individuals of M. lamarckii, 


divergence was conspecific 
which showed uncorrected divergence of 8.54%. 
Divergence could not be calculated for M. lyrata as it 
was represented by a single haplotype. Sequence 
variations between different species within the Meretrix 
ranged from about 1.83% to 14.81%. The uncorrected 
divergences of M. petechialis and M. lusoria with M. 
meretrix are 1.83%, 4.60%, which are lower the average 
uncorrected congeneric divergence of 11.08% of the 


genus Meretrix. (Tab. 3, Tab. 4). 

Among the 12 haplotypes, the third codon positions 
were invariant in only 83 out of the 191 sequenced 
samples in the Meretrix (Tab. 2). In contrast, more than 
four-fifths of the first codon positions (158 of 192) and 
about 95% of the second codon positions (184 of 191) 
were monomorphic. The number of parsimony 
informative sites at each of the codon positions for the 
whole data set, are also given in Tab. 2. There were 
approximately eleven times as many parsimony 
informative sites at the third codon positions (67) than 
the second codon positions (6), and about three times as 


many as at first positions (20). For the data set as a whole, 
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the ratios are 20: 6: 67 for first, second, and third 
positions, respectively—a ratio of about 3: 1: 11. To 
offset the preponderance of the third position variation in 
the reconstruction of phylogenetic trees, weighted 
analysis using a weighting of sites of 3: 11: 1 was 
included. 
2.2 Phylogenetic relationships 

Figs. 1-3 present the phylogenies recovered under 
parsimony and Bayesian analysis, respectively. The 
values of interior branch test for some of the nodes in the 
parsimony trees and the posterior probability values for 
some nodes in the Bayesian tree were low. However, the 
topologies of the trees were very similar in most clusters. 
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In this analysis, haplotypes of M. petechialis and M. 
lusoria are nested within M. meretrix forming clade I . 
The difference among all the trees is the position of M. 
lyrata. In parsimony tree (weighted), M. lyrata clustered 
with clade I, but is sister to the assemblage (clade I + M. 
lamarckii) in another parsimonious tree and Bayesian 
tree. 


3 Discussion 


3.1 Taxonomic position of M. petechialis and M. 
lusoria 

M. lusoria and M. petechialis were previously 

considered as independent species, different from M. 


Tab. 3 Intra specific nucleotide P distances for the Meretrix species 


Species 


Mean P distance 





1. M. meretrix 
2. M. petechialis 
3. M. lusoria 

4. M. lyrata 

5. M. lamarckii 


0.0325+0.0050 
0.0052+0.0030 
0.0697+0.0082 
n/c 
0.0854+0.0113 





Average intraspecific P distance in the genus = 0.045 


Tab. 4 Inter specific nucleotide P distances among the Meretrix species 





Species 1 2 3 4 5 
1. M. meretrix 

2. M. petechialis 0.0183 

3. M. lusoria 0.0460 0.0438 

4. M. lyrata 0.1407 0.1385 0.1440 

5. M. lamarckii 0.1448 0.1420 0.1414 0.1481 


Average interspecific P distance in the genus =0.1108 
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MP tree resulting from the analysis of the COI gene sequences of 12 haplotypes 


Numbers on nodes correspond to percentage bootstrap values for 1000 replicates. 
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Fig.2 MP tree resulting from analysis of the COI gene sequences of 12 haplotypes 
The tree is based on a weighting of codon positions 3:11:1 for first, second, and third codon positions, respectively. Numbers on nodes 


correspond to percentage bootstrap values for 1000 replicates. 
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Fig.3 Bayesian tree resulting from the analysis of the COI gene sequences of 12 haplotypes 


Numbers on nodes correspond to the value of posterior values. 


meretrix, in the genus Meretrix (Jukes-Browne, 1914; 
Habe, 1997; Zhuang, 2001). It was shown that the two 
species differed from M. meretrix by the following shell 
characteristics: Anterior and ventral margins, posterior 
end of shell. But the results from COI gene sequences 


analysis strongly disagreed with their viewpoints. 
According to the phylogenetic trees obtained by this 
study, M. petechialis and M. lusoria nested within M. 
meretrix, of them, M. meretrix has rounded posterior end, 
M. petechialis and M. lusoria are bluntly angled. 
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Therefore, the shape of the shell end is a homoplasious 
character, and species determination in the Mereterix 
using this character should be investigated further. The 
results obtained in this study do not support the present 
taxonomic status of M. lusoria and M. petechialis. This 
is in accordance with the ideas from allozyme analysis 
that M. lusoria and M. petechialis were the most closely 
related species within the genus Meretrix (Ayako et al, 
2008). This viewpoint also supported that M. Jusoria and 
M. meretrix belong to different geographic subspecies of 
one species (Pan et al, 2006). 

In our study, the COI gene fragment of M. lusoria 3 
is identical to M. petechialis 2, and M. petechialis 1 is 
identical to M. meretrix 2. Moreover, the unweighted 
divergences (1.83%, 4.60% and 4.38%) among these 
three 
interspecific divergence 


species are much lower than the average 
(11.08%). Prashad (1932) 
pointed out that M. meretrix is a species which 
experienced the greatest variation in the group of 
bivalves, and because of shades of shells and shell colors, 
it was wrongly divided into many species. Furthermore, 
Fischeer-Piette (1941) also suggested that M. petechialis 
and M. lusoria should be treated as the same species of 
M. meretrix. 

In summary, the views of taxonomic status of M. 
lusoria and M. petechialis are ambiguous based on 
morphology. In our study, the viewpoint that M. 
petechialis and M. lusoria should be treated as a junior 
synonym of M. meretrix was supported. 

3.2 Classification of the genus Meretrix 

The Meretrix is widely scattered through coastal 
areas of the Indian Ocean, Southeast Asia, China, Korean 
Peninsula and Japanese Archipelago. Because of the 
remarkable variation of shapes and patterns of shells, 
early researchers divided this genus into many species, 
such as M. meretrix, M. typical, M. petechialis, M. 
castanea, M. graphica, M. labiosa, M. inpudina, M. 
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